Partonic effects on the elliptic flow at relativistic heavy ion collisions 
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The elliptic flow in heavy ion collisions at RHIC is studied in a multiphase transport model. By 
converting the strings in the high energy density regions into partons, we find that the final elliptic 
flow is sensitive to the parton scattering cross section. To reproduce the large elliptic flow observed 
in Au+Au collisions at ^/s = 130A GeV requires a parton scattering cross section of about 6 mb. 
We also study the dependence of the elliptic flow on the particle multiplicity, transverse momentum, 
and particle mass. 
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I. INTRODUCTION 



Elliptic flow in heavy ion collisions measures the asym- 
metry of particle momentum distributions in the plane 
perpendicular to the beam direction. It results from the 
initial spatial asymmetry in non-central collisions [|l],^. 
Theoretical studies have shown that the elliptic flow is 
sensitive to the properties of the hot dense matter formed 
during the initial stage of heavy ion collisions (3|-|^ . 

Recently, elliptic flow has been measured at RHIC in 
Au-I-Au collisions at -^s = 130A GeV. A large elliptic 
flow of all charged particles near midrapidity was re- 
ported by the STAR collaboration Also, the trans- 
verse momentum dependence of the elliptic flow of all 
charged particles was measured by both the STAR |^] 
and the PHENIX |^ collaboration, while those of charged 
pions, kaons, protons and anti-protons were measured by 
the STAR collaboration (| jl^ . The experimental results 
show that the elliptic flow first increases with particle 
transverse momentum and then levels off. The depen- 
dences of elliptic flow on both the charged particle multi- 
plicity |^-§,|ll[ and the particle pseudorapidity have 
also been measured. 

To understand these experimental results, many theo- 
retical approaches have been used. These include semi- 
analytic models [p^, models with parton energy loss 




hydrodynamic models |^,|15|,^, transport mod- 
1|17| , and the hybrid model which combines hydrody- 
namic and hadronic transport models [ p8| . Among these 
studies, hydrodynamic models usually give the largest el- 
liptic flow and an almost linear increase in its value with 
the particle transverse momentum below 1.5 GeV/c. In 
the hybrid model of combining the hydrodynamic model 
with the RQMD transport model |l^] and choosing cer- 
tain effective equation of state, it is possible to obtain an 
elliptic flow that is comparable to the measured ones in 
heavy ion collisions at both SPS and RHIC energies [^ . 
In transport models including only the parton cascade, 
the elliptic flow has been shown to be sensitive to the 
parton scattering cross section, and a large value can be 
obtained with a large cross section ||2^,|2l| . On the other 
hand, transport models based on hadronic and/or string 



degrees of freedom in general give a smaller elliptic flow 
than that observed at RHIC. Including multi-Pomeron 
exchanges and hard gluon-gluon scatterings can, how- 
ever, yield a large elliptic flow [p7| . 

In this paper, we study the elliptic flow in Au+Au 
collisions at RHIC using a multiphase transport model 
(AMPT) that includes both initial partonic and final 
hadronic interactions as well as the transition between 
these two phases of matter p^-p^. In particular, we 
study the effect due to partons converted from the initial 
strings in high energy density regions. 



II. CONVERSION OF STRINGS TO PARTONS 

In most transport models for heavy ion collisions, such 
as the ART model [||, the RQMD model ||l|l, and the 
UrQMD model , initial primary collisions produce ei- 
ther hadrons or strings which later fragment to hadrons. 
In the HIJING model on the other hand, initial 

primary collisions also produce minijet partons which 
later enter into the string configurations and fragment 
to hadrons. These minijet partons are partons produced 
from initial hard collisions, i.e., from perturbative QCD 
processes involving a minimum transverse momentum 
transfer pg, which is chosen to be 2 GeV/c in the HI- 
JING model to reproduce the pp and pp data . 

In the AMPT model ||2|-|§, which takes initial condi- 
tions from the HIJING model, minijet partons first un- 
dergo scatterings before fragmenting into hadrons. Since 
the number of hard collisions in an A -I- ^ collision roughly 
scales as A^^^ and grows faster with colliding energy, 
while the number of strings roughly scales as A, minijets 
become more important in heavy ion collisions at higher 
energies. However, for central Au+Au collisions even at 
200 A GeV minijet partons account for only about 1/3 
of the total produced transverse energy, so the effect of 
parton scattering on the final particle multiplicities and 
spectra is quite small . 

The above picture of coexisting partons and strings 
during the initial stage of high energy heavy ion collisions 
is questionable when the energy density is much higher 
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than the critical density for the QCD phase transition. In 
this case, the strings are expected to melt into par tonic 
degrees of freedom. Both the transport model and 
the high density QCD approach [|9| predict that the ini- 
tial energy density of produced matter in central Au+Au 
collisions at RHIC is more than an order of magnitude 
higher than the critical energy density (~ 1 GeV/fm'^). 
Keeping strings in the high energy density region thus 
underestimates the partonic effects in these collisions. 

To model the string melting in high energy density 
regions to partons, we extend the AMPT model in the 
following way. After using the HIJING model (with jet 
quenching turned off) to produce the initial conditions, 
we first let strings fragment to hadrons using the LUND 
fragmentation [|0|j2|] buih in the PYTHIA routine |3|] 
and then convert these hadrons to partons according to 
their flavor and spin structure. In particular, a me- 
son is converted to a quark and an anti-quark, while a 
baryon is converted to three quarks, and an anti-baryon 
is converted to three anti-quarks, where quark masses are 
taken as the same as in the PYTHIA program e.g. 
ruu = 5.6 MeV, rud = 9.9 MeV, rus = 199 MeV. We 
further assume that quarks are produced isotropically in 
the rest frame of the hadron and start to interact only 
after a proper formation time Tq . In converting hadrons 
to partons, hadrons are not assigned a formation time as 
they are considered as an intermediate step in modeling 
the melting of strings to partons in an environment of 
high energy density. Based on the expectation that hard 
partons (e.g. those described by perturbative QCD) are 
produced early while soft partons are produced late in 
the process, we take Tq ~ l/Q for the parton proper 
formation time, where Q is a scale related to the par- 
ton transverse momenta. Since partons produced from 
a string through the same intermediate hadron should 
have the same formation time, the transverse mass of 
the hadron is thus used for Q in determining the parton 
formation time. 

The scatterings among these quarks are treated using 
the parton cascade ZPC with the following universal cross 
section |]32| : 

dap _ 97ra^ / 1 
dt 2 \ s ) [t- /i2)2 • ^ > 

In the above, the strong coupling constant as is taken to 
be 0.47, s and t represent the standard Mandelstam vari- 
ables for the two-parton elastic scattering process, and 
the effective screening mass /i depends on the tempera- 
ture and density of the partonic matter. In the present 
study, we shall take /i as a parameter to obtain the de- 
sired total cross section and the corresponding angular 
distribution. 

A parton can hadronize after it stops interacting, i.e., 
after it will no longer collide with other partons. We 
model the hadronization by combining the nearest two 
partons into a meson and three partons into a baryon (or 
an anti-baryon). As partons freeze out at different times 



and parton coalescence occurs at different times, the 
hadronization is treated locally. Since combinations of 
partons form a continuous invariant-mass spectrum but 
not a discreet one, it is in general impossible to conserve 
4-momentum when several partons are combined into a 
hadron. In our current treatment, we choose to conserve 
the 3-momentum and determine the hadron species ac- 
cording to the flavor and invariant mass of coalescing 
partons. E.g., if a u and a d quark coalesce, a 7r~ will 
be formed if the 2-quark invariant mass is closer to the 
7r~ mass, or a p~ will be formed if the 2-quark invariant 
mass is closer to the central value of the p mass |Q . All 
SU(3) mesons and baryons listed in the HIJING program 
|p7| are included for coalescence except 77',!]* and S*, 
which are not present in our hadronic transport model, 
and K^, states. The resulting hadrons are then given 
an additional formation time of tJ^ = 0.7 fm/c in their 
rest frame and imported to the ART hadronic transport 
model to take into account their rescatterings ||2^-|25|]. 

III. RESULTS ON ELLIPTIC FLOW 

A. Time evolution 

We have used the above extended AMPT model to 
study the elliptic flow in Au-|-Au collisions at RHIC en- 
ergies. The elliptic flow here is defined as 

^.^ _ hl-pl \ ^2) 

where the average is performed over all particles. For 
center-of-mass energy ^/s — 130A GeV and impact pa- 
rameter 6 = 8 fm, we show in Figs, ^a) and|l|(b) the time 
evolution of the elliptic flow of all partons, i.e., regard- 
less of their formation and freeze-out times, at midra- 
pidity with parton scattering cross sections Cp = 3 and 6 
mb, respectively. We note that partons still inside strings 
and partons which have frozen out are also included in 
Fig. 1. For both values of parton cross sections the ellip- 
tic flow is seen to develop mostly within the first 5 fm/c, 
and both the rate of increase and the final parton V2 are 
larger for the larger parton cross section. Also shown in 
these figures are the time evolutions of the parton trans- 
verse energy Et (scaled down by a factor of 2000), the 
average transverse energy per parton {Et) (scaled down 
by 2), and the second moment of the spatial asymmetry 
determined from the parton positions at their previous 
interaction points, i.e., S2 — {{x^ — y^^jix^ -I- y^))- It 
is seen that for both parton cross sections the absolute 
value of S2 decreases with time, and its final value is 
closer to zero in the case of cTp = 6 mb. The final satura- 
tion of parton vi in the partonic phase is due to the lack 
of scatterings in the limit of small and the vanishing 
spatial anisotropy in the limit of large Cp. We also see 
that the transverse energy at midrapidity decreases by ~ 
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30% within a couple of fm/c. Since the decrease of par- 
ton transverse energy is faster than that for the parton 
number at midrapidity, the average transverse energy per 
parton {Et) also decreases. The final saturation of {Et) 
at ~ 20% below the initial value reflects the equilibra- 
tion between the longitudinal and transverse momenta 
of partons HUl. 
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FIG. 1. Time evolutions of the parton elliptic flow, spa- 
tial anisotropy, transverse energy, and average transverse en- 
ergy per parton at midrapidity in Au-|-Au collisions at cen- 
ter-of-mass energy 130^4 GeV and impact parameter 6 = 8 fm 
for parton scattering cross sections = 3 mb (a) and 6 mb 
(b). 



In Fig. ^, we show the time evolution of V2 for active 
partons, i.e., partons that have not stopped scatterings, 
V2 for formed hadrons, the total V2 and S2 including both 
active partons and formed hadrons, and the number of 
active partons at midrapidity for collisions at an impact 
parameter 6 = 8 fm and a parton cross section cTp = 6 
mb. While the initial numbers of active partons are near 
zero because most partons are not yet formed and thus 
not yet active, it is seen that the number of active par- 
tons peaks within the first fm/c and decreases by a factor 
of 2 at about 4 fm/c due to the hadronization. The total 
V2 is thus dominated by active partons at the early stage 
and by formed hadrons during the later stage of heavy 
ion collisions. Although the active partons have a large 
elliptic flow at later times, their number is too small to 
affect the total elliptic flow. 
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FIG. 2. Time evolutions of elliptic flow, spatial anisotropy, 
and the number of active partons at midrapidity in Au-|-Au 
collisions at center-of-mass energy of 130A GeV and impact 
parameter fo = 8 fm for a parton cross section o-p = 6 mb. 



B. Impact parameter dependence 

We have also studied the impact parameter depen- 
dence of elliptic flow in Au-f Au collisions at center-of- 
mass energies of I30yl GeV. To compare with the cen- 
trality dependence in the STAR data Q, we first divide 
the impact parameter range < 6 < 13 fm into six bins 
with equal bin-size except the first bin, which is taken to 
be < & < 3 fm. The impact parameter dependence is 
then converted to the dependence on iVch/./Vmax by tak- 
ing Nch as the number of charged particles within the 
pseudorapidity range rj e (—0.75,0.75) and N^^^y, as its 
value at 6 = fm. Furthermore, we include only charged 
particles within rj G (—1.3,1.3) and transverse momen- 
tum range pt G (0.1, 2.1) GeV/c in evaluating the elliptic 
flow V2 in order to compare more directly with the STAR 
data on the centrality dependence. 

In Fig. H, we show the results for the elliptic flow 
of charged particles as a function of Nch/N^-ia-^ for the 
scenarios of the default AMPT (without string melting) 
and the extended AMPT (with string melting). The er- 
ror bars in our results represent only the statistical error 
in V2 but not that in Nch- Although all results show the 
qualitative features of the observed centrality dependence 
of V2 0], the shape and magnitude of V2 depend sensi- 
tively on the partonic dynamics JSGt . Without converting 
the initial strings into partons, the default AMPT model 
gives the smallest elliptic flow. Allowing the melting of 
strings to partons, the elliptic flow for a larger parton 
cross section not only is higher than that for a smaller 
partonic cross section, but also peaks at a lower value 
of Nch/Ntnax- Of the three parton cross sections, the re- 
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suits for dp = 6 mb appear to be more consistent with 
the observed centrahty dependence of the eUiptic flow. 
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FIG. 3. Impact parameter dependence of elliptic flow at 
130A GeV. The data from the STAR collaboration are 
shown by filled circles, while the theoretical results for differ- 
ent partonic dynamics are given by curves. 

We note that the results from the default AMPT model 
are insensitive to the parton scattering cross section, 
which is taken to be 3 mb in the scenario with strings. 
This is mainly due to the small fraction of energy that 
is carried by minijet partons and the lack of transverse 
collective motion of the strings. As a result, the elliptic 
flow is significantly reduced after minijet partons com- 
bine with strings and fragment to hadrons. However, 
with strings converting to partons, the initial energy orig- 
inally stored in the strings also contributes to the parton 
dynamics. This leads to a larger elliptic flow and its sen- 
sitivity to the parton cross section, and thus makes it 
possible to determine the strength of partonic interac- 
tions from the final elliptic flow. 



C. Transverse momentum dependence 

Our results for the dependence of charged particle el- 
liptic flow on the transverse momentum, i.e., the differen- 
tial elliptic flow V2{pt), are shown in Fig. ^. As observed 
in the experiment data, the differential elliptic flow first 
increases almost linearly with transverse momentum and 
then tends to level off at large transverse momenta. How- 
ever, both the slope of initial increase and the transverse 
momentum at which deviation from a linear dependence 
appears are affected by the parton dynamics. The result 
from the default AMPT model [p^-p^ (open triangles), 
which includes only minijets in the partonic phase, has 
the smallest V2{pt) at a given pt and shows a departure 



from the linear dependence also at the smallest transverse 
momentum. Including partons from string melting in the 
high energy density regions increases both the magnitude 
of 1)2 (pt) and the value of pt at which the linear depen- 
dence breaks down. 
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FIG. 4. Transverse momentum dependence of elliptic flow 
at 130A GeV. Circles are the STAR data for minimum-bias 
Au-I- Au collisions Q , and curves represent the minimum-bias 
results for charged particles within rj € (—1.3,1.3) from the 
AMPT model. 



The differential elliptic flow V2 [pt ) is expected to be 
different for different particles ||3^,|T^ . Fig. || shows our 
results for pions, kaons and nucleons from the extended 
AMPT model together with the STAR data for charged 
particles At low transverse momentum {pt < 1 

GeV/c), while the V2{pt) of pions and kaons increases al- 
most linearly with pt, that of nucleons shows a stronger 
dependence on pt- Furthermore, particles with smaller 
masses are seen to have higher values of V2{pt) at a given 
Pt- These features are qualitatively similar to those ob- 
served from the STAR data and also those obtained 
from the hydrodynamic models |l^,|l^. On the other 
hand, all particles seem to have similar values of elliptic 
flow, within the errors of our calculations, at a given pt 
above I GeV/c. 

To test the sensitivity of our results to different for- 
mation time parameters, we have calculated the elliptic 
flow at b = 8 fm for the following three cases for compar- 
ison: i) parton formation time Tq larger by a factor of 
2, ii) Tq smaller by a factor of 2, iii) hadron formation 
time set to O.Of fm/c instead of the default 0.7 fm/c. 
Compared to our default results, the relative changes of 
final charged particle V2 are found to be less than 10% 
for p^<2 GeV/c. 
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FIG. 5. Transverse momentum dependence of the elliptic 
flow for particles of different masses in the case of parton cross 
section ap — 6 mb in Au+Au collisions at ISOA GeV. Curves 
represent the minimum-bias results within ri G (—1.3, 1.3) 
from the extended AMPT model. 
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FIG. 6. Impact parameter dependence of elliptic flow 
at 200^1 GeV. Curves represent minimum-bias results for 
charged particles within rj G (—1.3,1.3) and pt G (0.1,2.1) 
GeV/c. 



D. Elliptic flow at ^ = 200 A GeV 

We have also studied the elliptic flow in Au-|-Au colli- 
sions at the maximum RHIC energy of 200A GeV using 
the extended AMPT model. In Figs. || and ^, we show 
the dependence of the charged particle elliptic flow on im- 
pact parameter and transverse momentum, respectively, 
together with the results for = 6 mb at y/s = 130 A 
GeV (curves with filled diamonds). Compared to Fig. ||, 
it is seen that for the same parton cross section the ellip- 
tic flow in less central events (A"ch/A^max < 0.6) increases 
only slightly (no more than 0.01 in magnitude) with the 
center-of-mass energy. As a result, the sensitivity of the 
centrality dependence of elliptic flow to the parton scat- 
tering cross section in heavy ion collisions at ^/s — 200A 
GeV is similar to that at ^/s = 130A GeV. 

For the differential elliptic flow. Fig. |^ shows that for 
transverse momenta below about 1 GeV/c, it shows even 
less change with the center-of-mass energy than the im- 
pact parameter dependence of the elliptic flow. As to the 
dependence on the parton cross section, the differential 
elliptic flow seems to show a larger sensitivity at higher 
transverse momenta. However, to study this quantita- 
tively requires much better statistics than we have ob- 
tained so far. Also, for elliptic flow at high pt, the effect 
of energy loss due to inelastic partonic processes, which 
have not been included in the AMPT model, becomes 
important jl^jlj]. Furthermore, the parton coalescence 
model may be less suitable than the independent par- 
ton fragmentation model for modeling the hadronization 
dynamics. 
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FIG. 7. Transverse momentum pt dependence of elliptic 
flow at 200 j4 GeV. Curves represent minimum-bias results 
for charged particles within rj G (—1.3, 1.3). 



IV. CONCLUSIONS AND DISCUSSIONS 



In this paper, we have studied the elliptic flow in 
heavy ion collisions at RHIC using a multiphase trans- 
port model that includes both partonic and hadronic 
scatterings. To take into account the effects of string 
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melting due to initial high energy density, we have intro- 
duced a schematic model to convert the strings produced 
in soft interactions into partons. We find that the magni- 
tude of the elliptic flow is sensitive to the scattering cross 
sections of these partons. To reproduce the large elliptic 
flow observed in Au-|-Au collisions at 130^ GeV at RHIC 
requires a parton scattering cross section of roughly 6 
mb. Similar to the findings from hydrodynamic models, 
the differential elliptic flow V2{pt) for charged particles is 
found to increase almost linearly with pt at low transverse 
momentum. At high transverse momentum V2{pt) devi- 
ates from a linear dependence and becomes more flat. 
Also, heavier particles have smaller V2{pt) than lighter 
particles at a given low pt, while they seem to have sim- 
ilar values of elliptic flow at higher pt. We further find 
that the increase of the elliptic flow from — 130A to 
200A GeV is quite modest in the centrality dependence 
and is even less in the dependence on transverse momenta 
below 2 GeV/c. 

In the present study, we have adopted a simple ap- 
proach in converting strings to partons and in the 
hadronization from partons to hadrons. To treat this 
more consistently, we need to consider the local par- 
ton and string densities and determine when strings are 
converted to partons and partons are combined to from 
hadrons. In a probably more systematic approach, the 
initial energy may be separated into two parts with the 
hard part corresponding to minijet gluons and quarks, 
and the soft part corresponding to color fields. A par- 
ton cascade including a color mean field can then be 
used to describe the subsequent dynamics. Furthermore, 
both inelastic parton scatterings and particle subdivi- 
sion [ p8|j2l]j39[ | have been neglected in the parton cascade. 
These may affect the elliptic flow and need to be studied 
in the future. 
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Li, T.W. Ludlam, M. Murray, J.Y. OUitrault, Subrata 
Pal, D.H. Rischke, N. Xu and B. Zhang. This work 
was supported by the U.S. National Science Founda- 
tion under Grant Nos. PHY-9870038 and PHY-0098805, 
the Welch Foundation under Grant No. A-1358, and 
the Texas Advanced Research Program under Grant No. 
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